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ABSTRACT 

Unveiling the nature of cosmic dark matter is an urgent issue in cosmology. Here 
we make use of a strategy based on the search for the imprints left on the CMB 
temperature and polarization spectra by the energy deposition due to annihilations of 
the most promising dark matter candidate, a stable WIMP of mass = 1 — 20 GeV. 
A major improvement with respect to previous similar studies is a detailed treatment 
of the annihilation cascade and its energy deposition in the cosmic gas. This is vital 
as this quantity is degenerate with the annihilation cross-section (cru). The strongest 
constraints arc obtained from Monte Carlo Markov Chains analysis of the combined 
WMAP7 and SPT datasets up to £max = 3100. If annihilation occurs via the e+ — e~ 
channel, a light WIMP can be excluded at 2-a c.l. as a viable DM candidate in the 
above mass range. However, if annihilation occurs via ji^ — or — t~ channels 
instead we find that WIMPs with > 5 GeV might represent a viable cosmological 
DM candidate. 

We compare the results obtained in the present work with those obtained adopt- 
ing an analytical simplified model for the energy deposition process widely used in 
literature, and we found that realistic energy deposition descriptions can influence the 
resulting constrains up to 60%. 

Key words: dark matter 



1 INTRODUCTION 

According to the widely accepted Lambda Cold Dark Mat- 
ter (ACDM) cosmology, the Universe is mostly made of dark 
components, i.e. Dark Energy (75% of the mass-energy bud- 
get) and Dark Matter (DM, 20%); these components largely 
dominate over baryons (Komatsu ct al. 2011). The situation 
is then rather unsatisfactory as the nature of the dark com- 
ponents is far from being established and it stands as one of 
the most crucial issues in cosmology. 

The most promising DM interpretation is in terms of 
a thermal relic density of stable weakly interacting massive 
particles (WIMPs). An appealing feature of such a scenario 
is that the annihilation cross-sections predicted by the elec- 
troweak scale automatically provide the right DM density 
after freeze-out (Bertone & Silk 2010). This argument ap- 
plies equally well to particles with 1 — 20 GeV masses as to 
those with masses more traditionally associated with super- 
symmetric neutralinos (m^ ~ 40 — 1000 GeV). 

In the recent years, evidences have been accumulating 
in favor of DM in the form of ~ 10 GeV WIMPs. In fact, 
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a relatively light DM particle with an annihilation cross- 
section consistent with that predicted for a simple ther- 
mal relic {{av)T ~ 10~^® cm^ s~^) and a distribution in 
the halo of the Milky Way consistent with that predicted 
from simulations could accommodate the indirect detection 
of gamma rays from the Galactic Center, the synchrotron 
emission from the Milky Way radio filaments and the dif- 
fuse synchrotron emission from the Inner Galaxy (the so- 
called "WMAP Haze") (Hooper, Finkbeiner & Dobler 2007; 
Hooper & Linden 2011a, b; Dobler et al. 2010; Finkbeiner 
2004). 

At the same time it would be compatible with 
claims of low-energy signals from DM direct detection 
experiments as DAMA/LIBRA, CoGeNT, and CRESST- 
II. In particular, the striking detection of annual mod- 
ulation observed by DAMA/LIBRA (now supported by 
CoGeNT) is hardly consistent with all known stan- 
dard backgrounds. Notice, however, that (a) other ex- 
periments, as XenonlOO, have not confirmed the result 
of the direct detections, and (b) indirect detection fea- 
tures might have alternative astrophysical explanations 
(Bernabei et al. 2008; Bicrmann et al. 2010; Akcrib ct al. 
2010; CDMS II Collaboration et al. 2010; Bernabei et al. 
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2010; Crocker & Aharonian 2011; Aalseth et al. 2011a,b; 
Aprile et al. 2011; XENONIOO Collaboration et al. 2012; 
Guo & Mathews 2012). 

A phenomenological model of light DM particle able to 
accommodate the collection of indirect and direct observa- 
tions should require that DM annihilates primarily into lep- 
tons with a cross-section close to {av)T- Moreover, approx- 
imately 20% of annihilations must also proceed to hadronic 
final states in order to yield a spin-independent, elastic scat- 
tering cross-section (~ 10"''^ cm'^) with nucleons compatible 
with the direct detection (see Hooper 2012, for a detailed re- 
view) . 

The light DM hypothesis implies a larger cosmic num- 
ber density of such particles (wdm oc Q,-DMh^ /iti-dm)', in ad- 
dition, the annihilation rate (cx oc (1 + 2)") increases 
dramatically at early cosmic times. These two facts imply 
that the annihilation energy deposition might profoundly af- 
fect the thermal and ionization history of the Intergalactic 
Medium (IGM)^ prior to reionization. In turn, this modi- 
fied evolution with respect to the standard recombination 
scenario can in principle leave detectable signatures on the 
CMB spectrum^. Determining the amplitude of this effect 
is the chief goal of the present study. 

The effects of the DM annihilation around the red- 
shift of the last scattering surface (LSS) have been discussed 
in Slatyor, Padmanabhan & Finkbcincr (2009) and are only 
briefly summarized here. The extra free-electrons resulting 
from the DM energy cascade scatter CMB photons, thus 
thickening the LSS and in principle shifting the position of 
the peaks in the temperature-temperature (TT) power spec- 
trum. In practice, reasonable electron density excesses yield 
corrections to the positions of the peaks that can be safely ig- 
nored here. More importantly, oscillations on scales smaller 
than the LLS width are damped in the TT and EE spec- 
tra in a manner inversely proportional to their wavelength. 
Such DM annihilation effects on the TT spectrum are degen- 
erate with variations of the slope (tIs) and amplitude (As) 
of the primordial power spectrum, and, to a lesser extent, 
with the baryon (Qbh^) and DM (r^DM^i^) density parame- 
ter. Polarization spectra are generated via Thomson scatter- 
ing of the local quadrupole in the temperature distribution. 
As the broadening of the last scattering surface increases 
the intensity of the quadrupole moment, the EE spectrum 
is enhanced on large scales. Furthermore, it can be shown 
(i.e. Padmanabhan & Finkbcincr 2005) that the quadrupole 
is dominated by the free-streaming from the dipole pertur- 
bation that is 7r/2 out of phase of the monopole. A thicker 
LSS boosts the fractional contribution from the monopole, 
thus slightly shifting the peaks of the EE and TE spectra. 

A key aspect of these calculations is that only a 
fraction of the released energy is finally deposited into 
the IGM in the form of heating and H/He ionizations. 
However, earlier studies (Padmanabhan & Finkbeiner 2005; 
MapelU, Ferrara & Picrpaoli 2006; Galli ct al. 2009) have 
used a simplified description of such processes, based 

^ Strictly speaking, the term "Intergalactic Medium" is ill- 
defined before the epoch of galaxy formation. Nevertheless, fol- 
lowing common practice, we will use it anyway. 
^ Leptons and photons deposit their energy very efficiently in the 
IGM; hence leptonic channels are expected to leave a stronger 
CMB signature 



on the hypothesis that a redshift-independent fraction 
of the DM rest-mass energy is absorbed by the IGM. 
More recently (Slatyer, Padmanabhan & Finkbcincr 2009; 
Galli ct al. 2011; Hiitsi et al. 2011) have reassessed the 
energy deposition problem including various energy-loss 
mechanisms in a more realistic way. This approach, based 
on semi-analytical solutions lacks an implementation of 
low-energy atomic processes that determine the actual 
absorption channel (e.g. heating, ionization, excitations) 
and because of this they have to rely on the results 
of Chen & Kamionkowski (2004). To fill this gap here 
we build upon our previous work (Valdes, Evoli & Ferrara 
2010) in which we developed the MEDEA (Monte Carlo Energy 
Deposition Analysis) code which includes Bremsstrahlung 
and inverse Compton processes, along with H/He coUi- 
sional ionizations and excitations, and electron-electron col- 
lisions. MEDEA enables us to compute the energy partition 
into heating, excitations and ionizations as a function of 
the primary initial energy, the gas ionization fraction and 
the redshift. MEDEA has been recently improved (Evoli et al. 
2012) to include the energy cascade from particles gen- 
erated by primary leptons/photons using the most up-to- 
date cross-sections and extending the validity of the model 
to unprecedented high (~ TeV) energies (see ShuU 1979; 
ShuU & van Steenberg 1985; Furlanetto & Stoever 2010). In 
addition, arbitrary initial energy distribution of electrons, 
positrons and photons can be assigned. These improvements 
make MEDEA suitable for studying the IGM energy deposition 
for some of the most popular DM candidates (Evoli et al. 
2012). With this greatly improved physical description we 
aim at computing the signatures left in the CMB spectrum 
by annihilating light DM. 



2 METHOD 

In this Section we compute the energy input of DM 
annihilations in the IGM. This approach is similar in spirit 
to a number of recent works (Padmanabhan & Finkbeiner 
2005; GaUi et al. 2009; Hiitsi, Hektor & Raidal 2009; 
Slatyer, Padmanabhan & Finkbeiner 2009; Galli et al. 
2011; Hiitsi et al. 2011); however, we improve upon them 
by a more accurate description of the energy deposition 
channels. 



2.1 Modified ionization history 

For the reasons given in the Introduction, we concentrate 
on light DM candidates that annihilate mainly in leptonic 
channels. In Fig. 1 we show the annihilation spectra of a 
10 GeV DM particle for the different annihilation channels, 
computed using the public code DarkSUSY. The muonic and 
tauonic channels produce a leptonic pair whose prompt an- 
nihilation gives rise to an energy spectrum of primary elec- 
trons or positrons with kinetic energy from 10 GeV down to 
few tens of MeV; annihilation in the electron channel pro- 
duces a electron/positron pair in which both the two pri- 
mary leptons have a kinetic energy which is the mass of the 
annihilation particle. 

The total energy density input from DM annihilations 
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considered here as it depends strongly on the DM model 
chosen (van den Aarssen, Bringmann & Goedecke 2012). 

To derive the DM-modified cosmic ionization/thermal 
history, we need to include the above heating (and corre- 
sponding ionization) rate into the relevant detailed balance 
equations. To this aim we have modified that publicly avail- 
able code' RECFAST (Seager, Sasselov & Scott 1999), part of 
the CAMB (Lewis, Challinor & Lasenby 2000) code, by adding 
the following terms: 



Figure 1. Energy spectrum of electrons or positrons from the 
annihilation of a 10 GeV mass WIMP into e+e~, and 
T+r^ channels. 
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where pc — 3Hq /SttG is the critical density of the universe 
today, fioM is the DM density contribution to the critical 
density, ttidm is the mass of the DM particle and (av) is the 
thermally averaged product of the cross-section and relative 
velocity of the annihilating DM particles. Moreover we have 
defined B = (av) /3 x lO"^** cm^ s"\ Notice that eq. 1 is 
valid only for DM Majorana particles. 

In the light of the earlier works of Cirelli, locco & Panci 
(2009) and Hiitsi et al. (2011) we neglect the role of struc- 
ture formation in the calculation of the energy deposition. In 
fact, halos with density higher than the background could in 
principle boost the average annihilation rate, however their 
formation starts at a relative low-redshift {z 100) when 
the ionization rate due to DM annihilation is already negli- 
gible. 

By introducing the mean number density of hydrogen 
nuclei nn ~ 1.9 x 10~'^(1 ~f zf 



It has been pointed out that (av) could be somewhat 
boosted by the Sommerfeld effect (e.g. Galli et al. 2009; 
Slatyer, Padmanabhan & Finkbeiner 2009). Although it is 
easy to implement this process in this scheme we have not 
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where /ion,H, /ion, Ho (/ii) are energy deposition fractions into 
H or He ionizations (heating) including those induced by 
Lya photons on atoms in the excited states. 

A key point to take from eqs. 4—6 is that energy de- 
position fractions are fully degenerate with the parameter 
we aim to constrain, i.e. {crv). To partly alleviate this diffi- 
culty, a possible strategy, first proposed by Ripamonti et al. 
2005, is to determine the lowest possible bound by assuming 
fh = /ion = 1. More often, constraints have been derived by 
using the Chen & Kamionkowski (2004) prescription for fi. 
Based on the results of ShuU & van Stccnbcrg (1985), these 
authors pointed out that when the gas is mostly neutral, en- 
ergy is evenly distributed among ionizations, excitation and 
heating; for a fully ionized medium, almost all of the energy 
goes instead into gas heating. A linear interpolation is used 
for intermediate ionization values: 
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This approximation is too crude to be used for high- 
precision predictions as clearly shown by the comparison 
with fully-fledged Monte Carlo simulations (Valdcs ct al. 
2007; Furlanetto & Stoever 2010). Moreover, for primary en- 
ergies ^ 1 MeV inverse Compton energy losses on the CMB 
becomes important and introduce a significant redshift de- 
pendence of the fractions. These processes have been care- 
fully modeled in Evoli ot al. (2012) and here we use their re- 
sults for fi. Note that the latter assume that photon energy 
deposition occurs locally, which is not true in general (see 
Slatyer, Padmanabhan & Finkbeiner 2009). In Appendix A 
we show that this approximation is very accurate in the 
energy range of interest here. The energy depositions calcu- 
lated with MEDEA are shown in Fig. 2 for different annihi- 
lation channels of a 10 GeV DM particle mass. Such curves 
show a dependence on the annihilation channel since differ- 
ent initial spectral distributions involve different energy loss 
mechanisms. For computational speed-up purposes we have 



Recently other similar codes (e.g. CosmoREC, HyREC) have 
improved the precision of the results implementing a more de- 
tailed description of the atomic structure (Hiitsi et al. 2011; 
Giesen et al. 2012). However, given the current precision of CMB 
experiments, these corrections do not affect our conclusions. 
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Figure 2. Fractional energy depositions into H and He ionization 
(solid lines) and heating (dotted lines) from DM annihilation of 
a 10 GeV DM candidate as function of redshift. Color code as in 
Fig. 1. The orange lines show the energy deposition fractions ob- 
tained from the simplified deposition model described in Sec. 3.1 
for annihilation into muons. In this plot we assume the standard 
IGM ionization evolution based on a WMAP7 cosmology. 

derived handy fitting formulae, given in Appendix B, to the 
MEDEA numerical results. 

2.2 MCMC analysis 

To obtain a constraint on the annihilation cross-section of 
light DM candidate, we have performed a Monte Carlo 
Markov Chains (MCMC) analysis using publicly available 
CosmoMC package (Lewis & Bridle 2002). We consider here 
a flat ACDM model with the canonical six parameters plus 
an additional seventh one, {(jv). Therefore the theoretical 
model we adopt is described by the following set of param- 
eters: 

{u}b,U)'Dyi,Os,T,ns,\og[lQ^° As], {(Tv)}, (9) 

where = Qth^, uj-dm = ^lomh^ are the baryons and cold 
dark matter density parameters, 6s is the ratio between the 
sound horizon and the angular diameter distance at decou- 
pling, r is the optical depth, Us is the scalar spectral index. 
As is the amplitude of the primordial spectrum. The flat 
priors assumed for these parameters are shown in Tab. 1. 

Our basic data set is the seven-year WMAP tem- 
perature and polarization data (Komatsu ct al. 2011; 
Larson et al. 2011). We consider purely adiabatic initial con- 
ditions and we impose spatial flatness. We also fixed the 
primordial fractional abundance of Helium to the standard 
observed nominal value of Yne ~ 0.24. We refer to this ba- 
sic dataset as "WMAP7" . For each case we run 5 chains; 
convergence diagnostic tests are performed using the Gel- 
man and Rubin "variance of chain mean/mean of chain vari- 
ances" R-1 statistics. We consider the chains to be converged 
only if i? - 1 < 0.03. The 68% and 95% confidence level 
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Table 1. Adopted flat priors for the cosmological parameters. 



(c.l.) one- and two-dimensional constraints are obtained af- 
ter marginalization over the remaining "nuisance" parame- 
ters. We have tested that varying Ho instead of 6s, as sug- 
gested in Galli ct al. (2009), our results are affected by less 
than 5%. 

In addition to the WMAP7 dataset we also consider the 
case ^^CMB ALL+SPT" . In this larger dataset we include, 
in addition to the WMAP data, the CMB temperature and 
polarization data from QUaD (Brown ct al. 2009), and the 
recent SPT (Keisler et al. 2011) data. The inclusion of the 
QUaD experiment (a) enlarges the multipole range consid- 
ered for the temperature, allowing to probe the small-scale 
region 500 ^ ^ 2500, and (b) adds information on the E- 
and B-mode polarization. Moreover, the SPT experiment 
pushes the dynamic range of CMB observations to larger 
multipoles with the respect of WMAP7, measuring with a 
better accuracy the damping tail of the CMB angular power 
spectrum. We consider data up to £ = 3100. For the SPT 
experiment it is necessary to account for foreground contri- 
butions by adding three extra parameters representing the 
amplitude of the SZ, Asz, clustering, Ac, and shot-noise, 
Ap, signal from point sources. We used for each foreground 
component the proper template provided by (Keisler ct al. 
2011). When deriving our constraints we marginalize over 
these three nuisance parameters. To compute the likelihood 
of the data we have properly modified the CosmoMC pack- 
age in order to make use of the routine supplied by the 
WMAP team for the WMAP7 dataset, publicly available at 
the LAMBDA website* , and of the likelihood code provided 
by the SPT team (Keisler et al. 2011) for the SPT experi- 
ment. 

As we already discussed in the introduction the inclu- 
sion of small scale CMB measurements can greatly help in 
breaking the degeneracy with the other cosmological param- 
eters, and in particular with Us, thus improving the con- 
straints on the DM sector parameters. Moreover, the ad- 
dition of the SPT data to the WMAP data improves the 
constraints on the the ratio of the sound horizon to the an- 
gular diameter distance parameter ds by nearly a factor of 
two (Keisler et al. 2011), thus narrowing the allowed range 
of the other parameters. 

We adopt the standard parametrization for the reion- 
ization, considered as an instantaneous process occurring 
at some redshift Zr, with Zr < 32. Such choice leads to a 
one-to-one relation between Zr and the adopted e.s. opti- 
cal depth T. As a caveat, we note that Pandolfi et al. (2011) 

4 http://lambda.gsfc.nasa.gov/ 
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Figure 3. DM-modified ionization (top) and thermal (bottom) 
histories for a 10 GeV WIMP annihilating into muons. Values 
of the annihilation cross-section correspond to different curves as 
shown by the legend. The black solid line represents the case case 
without DM annihilations. 
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Figure 4. Angular power spectrum of CMB temperature fluc- 
tuations: standard case without DM annihilations (black line), 
considering the /^+ — /^+ a 10 GeV WIMP with (av) = 3 X 
10"^*^ cm'^ s~^ (red line). The points with errorbars show the 
7-year measurements of the WMAP satellite (black) and the SPT 
data (blue). 



showed that a more realistic reionization modeling might af- 
fect the cosmological parameters that are more degenerate 
with the DM annihilation cross-section, thus introducing an 
additional source of uncertainty. 



3 RESULTS 

DM-modified ionization and thermal histories for a 10 GeV 
WIMP annihilating into muons on top of a ACDM model 
are shown in Fig. 3; the corresponding TT (EE) spectrum 
is shown in Fig. 4 (Fig. 5). Qualitatively similar conclusions 
can be drawn for the other channels. The energy released in 
form of electrons and positrons from the annihilation of DM 
particles delays and quenches the recombination processes, 
thus resulting in a freeze-out relic electron fraction a factor 
of a few larger, depending on the value of (ctv). For the same 
reason, the temperature drop with time is less pronounced. 
As a consequence of the higher ionization rate, the CMB 
normalization value is smaller. 

We follow the procedure described in Sec. 2.2 to get 
constraints on the cosmological parameters in (9) and we 
compare them with those obtained by the WMAP team 
from their 7-yrs data. We present in Fig. 6 the 2-sigma 
c.l. constraints on the DM annihilation cross-section {(jv) 
as function of the DM mass. Differently to (e.g. Galli et al. 
2009) our results cannot be given as a single number due 
to the mass dependence of the energy deposition fractions 
(see Sec. 2.1). A detailed comparison with their results will 
be given in Section 3.1. The main conclusion is that only 
DM candidates lighter than 10 GeV annihilating via the 
e"*" — e~ channel can be excluded as a dominant component 
of the DM energy-density. The constraints are stronger, as 
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Figure 5. As in Fig. 4 for the polarization fluctuations. The points 
with error bars show the measurements of the QUaD experiment. 



expected, if we include in the present analysis the recent 
SPT dataset with ^max = 3100^ and the polarization data. 
In this case the electron channel is excluded in the entire 
mass range (up to 20 GeV), where the other two channels 
can be excluded for masses ^ 5 GeV. 



We have verified that our results are insensitive for choice of 
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We have verified that the stronger constraints come 
mainly from the SPT data inclusion, since the polariza- 
tion data alone improve the constraints by < 3%. Cur- 
rently polarization data alone are not of sufficient quality 
to robustly constrain DM parameters. Future experiments, 
specifically devoted to measure polarization at smaller scales 
like Planck (Tauber et al. 2010), PolarBear (Anthony 2012) 
and CMBPol (Zaldarriaga et al. 2008) are expected to sig- 
nificantly improve the situation. 

The CMB constraints we find are weaker than the 
constraints obtained by the FERMI experiment using the 
signal in the diffuse isotropic gamma emission from the 
Galax;y (Abdo et al. 2010) and from a combined analy- 
sis of the Milky Way satellites Ackermann ct al. (2011); 
Baushev, Federici & Pohl (2012); Cholis & Salucci (2012). 
Comparing the 10 GeV case of annihilation channel in 
muons and taus, the inferred maximum cross-section from 
FERMI falls below the thermal value. However, in their anal- 
ysis the rather uncertain distribution of DM in galaxies must 
be specified, while the present approach is free from any such 
hypothesis. 

In Table 3 we report the 68% c.l. constraints on the 
cosmological parameters for the 10 GeV muon annihilation 
channel fot the WMAP7, CMB ALL+SPT cases, and the 
WMAP7 alone dataset, ie. a minimal ACDM model with- 
out annihilating DM ("WMAP7 (Standard)"). The one- 
dimensional posterior probability for flbh^, floMh^ and Us, 
for the three dataset cases considered is also shown in 
Fig 7. The strongest shift occurs for the baryon density 
Qf,h^ which in the minimal, 6-parameter, standard case is 
Qbh^=0.0226 ± 0.0006, whereas, after the inclusion of the 
annihilating DM, becomes r26fe^=0.0224±0.0006 in the case 
of WMAP7, and ^6/1^=0.0217 ± 0.0004 in the case of CMB 
ALL+SPT. This lower baryon density required results from 
the increased number of electrons produced DM annihila- 
tions; the two factors combine to give the same optical fi- 
nal depth needed to match the CMB data. The constraints 
on the dark matter density are only barely affected by the 
introduction of the DM annihilation, while instead the con- 
straints on the scalar spectral index of primordial pertur- 
bations are shifted to higher values. Similarly to the case 
of ftbh^, but in the opposite direction, the extra energy in- 
jected by the DM annihilation leads to a damping of the tail 
of CMB power spectrum, so that Us has to be increased in 
order to compensate for this effect and still provide a good 
fit to the data. Notice that in the case of WMAP7, the intro- 
duction of DM annihilation makes the Harrison-Zel'dovich 
value for the scalar spectral index Us — 1 compatible with 
the data within two standard deviations, while instead when 
also the SPT dataset is added the scale invariant power spec- 
trum is again ruled out by the data. 



WMAP7 



3.1 Simplified energy deposition model 

As we have stressed already, using a correct description of 
the energy deposition fractions is crucial to derive reliable 
DM constraints. Here we intend to quantify this statement 
by comparing our results with the constraints obtained using 
an approximated energy deposition model. 
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Figure 8. Relative differences [(MEDEA-simplified)/MEDEA] 
between tlic constraints obtained with a simple energy deposition 
model as described in Sec. 3.1 and using the energy deposition 
fractions obtained with the MEDEA code. Color code as in Fig. 1 



This is summarized by the following expressions: 
^ 1- Xe /=,1 + 2/hc - X 
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where x = xh + /Hea;He is a convenient variable to be used 
in RECFAST and 



Ch = Ch + (1 - Ch) 



-Bion.H 
Ea.U 



(13) 



(a similar expression is valid for the He) where 
Ch and Chg are the Peebles factors as given 
in Wong, Moss & Scott (2008). As in Galli et al. (2011) 
we have multiplied these formulas for the /abs(z) given 
by Slatyer, Padmanabhan & Finkbcincr (2009) for the DM 
annihilation in electrons or muons at 1 GeV. 

In Fig. 2 we show the corresponding energy depositions 
as function of redshift for the muon channel and we compare 
with what obtained from the Monte Carlo simulations. It is 
evident that this simplified approach over predict the energy 
deposition for almost the entire redshift range. 

We have verified that using the analytic expression in 
eq.s 10-12 the derived constraints at 1 GeV are coincident 
with the results reported in Tab. II by Galli et al. (2011) 
either for the muon and the electron channel. 

In Fig. 8 we show the relative differences between our re- 
sults and the results obtained adopting the simplified model. 
We compare the case in which only WMAP7 data are used. 
In the range mDM = 1 — 20 GeV the differences can be 
quoted between 10% and 30% for the electron channel, and 
between 20% and 60% for the muon channel. The con- 
straints we get always tend to be weaker than those given 
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Figure 6. Constraint plot on the maximum cross-section for different DM candidates based on WMAP7 and CMB ALL-j-SPT dataset 
(the region excluded for the tau annihilation channel is indicated in green, while the additional region excluded for the muon (electron) 
annihilation channel is indicated in red (blue)). 
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Figure 7. One dimensional posterior probability distribution of Q,i,h? , finivi^^ and ris parameters in the case of WMAP7 dataset and 
no DM annihilation (dashed), WMAP7 case (dotted) and CMB ALL + SPT case (dotted) 

Table 2. 95% c.l. upper limit constraints on the DM annihilation cross-section {crv) [10~^® cm^/s] in different cases of mass, annihilation 
channel and dataset considered. 



Mass 




WMAP7 




CMB ALL+SPT 




e+e 




T+r 


e+e 




T+r 


2 GeV 


^0.554 


^2.19 




sJO.351 


^1.02 




5 GeV 


s;i.29 


s;5.07 


^6.91 


^0.811 


^2.55 


^3.80 


10 GeV 


<2.66 


<8.79 


^21.2 


^1.68 


^5.14 


^11.8 


20 GeV 


^6.32 


^16.2 


^40.0 


^3.60 


^10.1 


^26.7 
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Parameter WMAP7 (Standard) 



WMAP7 



CMB ALL+SPT 



ris 

In (10^0 A, 



0.0226 ± 0.0006 
0.1109± 0.0056 
1.0388± 0.0027 
0.0884± 0.0152 
0.9635± 0.0142 
3.1904± 0.0457 



0.0224±0.0006 
0.1105± 0.0054 
1.0384± 0.0027 
0.0876 ± 0.0150 
0.9716 ± 0.01479 
3.1868±0.0460 



0.0217± 0.0004 
0.1099 ± 0.0049 
1.0397±0.0015 
0.0831± 0.0138 
0.9579± 0.0113 
3.2122 ± 0.0446 



Table 3. 68% c.l. constraints of the background cosmology parameters in the case of WMAP7 datasct with no DM annihilation (WMAP7 
(standard)), compared with the WMAP7 case with DM annihilation, and the CMB ALL+SPT case. 



by Gain et al. (2011): the difference originates from the in- 
clusion of the low-energy processes inducing a net energy loss 
(i.e energy not going into heating, ionization or excitation). 
As explained in the Introduction, decreasing the energy de- 
position fractions makes the constraints weaker. 



4 SUMMARY AND DISCUSSION 

We have investigated the imprints left on the CMB temper- 
ature and polarization spectra by the energy deposition due 
to annihilations of one of the most promising DM candi- 
date, a stable WIMP of mass = 1 — 20 GeV annihilating 
into leptons. A major improvement with respect to previous 
similar studies is a detailed treatment of the annihilation 
cascade and its energy deposition in the cosmic gas. This is 
vital as this quantity is degenerate with (av). 

We performed a Monte Carlo Markov Chain analysis us- 
ing a modified version of the CosmoMC code and CMB data 
from the WMAP, QUaD and SPT experiments. By further 
marginalizing over the cosmological parameters of the back- 
ground cosmology we obtain the constraints on the annihi- 
lation cross-section for each annihilation channel. 

The strongest constraints are obtained by combining all 
the available datasets up to ima,x = 3100. If annihilation oc- 
curs via the e+ — channel, a light WIMP can be excluded 
as a viable DM candidate in the above mass range. How- 
ever, if annihilation occurs via jj.'^ — /i^ or r"*" — t~ channels 
instead we find that WIMPs with rrix > 5 GeV cannot be 
ruled out at 2-a c.l. to provide the cosmologically required 
DM content. 

We have compared our results with the constraints ob- 
tained by assuming a simplified energy deposition model, 
as the one profusely used in the recent literature, and we 
found that realistic energy deposition descriptions can infiu- 
ence the resulting constrains up to 60%. A different analysis 
of the energy deposition uncertainties has been performed 
in Gain ct al. (2012). 

We expect that this approach will be relevant in partic- 
ular with the upcoming Planck data, with its better sensitiv- 
ity and which allow a better constraining of this additional 
sources of ionization. 
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APPENDIX A: LOCAL DEPOSITION 

We have assumed in this paper that photon energy de- 
position occurs locally, which is not true in general (see 
Slatyer, Padmanabhan & Finkbeiner 2009). In the following 
we show that this approximation is very accurate in the en- 
ergy range of interest here. CMB photons gain energy as 
they are Inverse Compton scattered by energetic leptons. 
At each scattering event, a CMB photon with mean energy 
i?7,CMB will be upscattered to an energy equal to: 

where 7 is the Lorentz factor for the lepton. At 
epochs in which energy deposition is important 
{z ^ 1000) such upscattered photons are subse- 
quently mainly downgraded by Compton scattering 
with thermal electrons (Chen & Kamionkowski 2004; 
Slatyer, Padmanabhan & Finkbeiner 2009). 

To estimate the efficiency of this mechanism, we 
compare in Fig. Al the Compton cooling time, t~^^^ = 
(dlnE/dt) — carrieil + z)^eg{e), with the Hubble time. 



20 




z 



Figure Al. Ratio between Hubble time and Compton cooling 
time scales as function of the redshift. 

tn = H~^{z) where e = E^/m^c? and (7(e) is the classical 
Klein-Nishina cross-section. It is evident that the local de- 
position assumption, requiring tn ^ ^cooi, can be safely ap- 
plied in the redshift range of interest here (600 ^ z ^ 1000) 



APPENDIX B: FITTING FORMULAE 

Below are the fitting formulae to the numerically derived 
energy depositions of electrons and positrons in the various 



channels. 

fH{x.,z) = 10^(^)(l-C(2)(l-a;f(^))) (Bl) 

hon,H{^e,z) = 10^^^^ ( 1 - X^^^^ ) ^ ^ (B2) 

fior.M^.,z) = 10^('-)(l-x^'^')^'^) (B3) 
where: 

A{z) = Ao-F Ailogio2H-42(logio2)' (B4) 

B{z) = Bo + Bilogioz + B2(logio2)^ (B5) 

C{z) = Co + Cilogio2-f C2(logio^)' (B6) 



The values of the parameters are given in Tab. Bl. 

Moreover we provide an updated version with respect 
to Evoli et al. (2012) for the energy deposition fractions that 
can be used for energies below the IC threshold: _Eth = ((l-f 
z)/21)-i/2 MeV: 

fh{x^) = a{l~c{l^xl)) (B7) 
/ion(a-e) = a{l-xlY (B8) 
where the values of the parameters are given in Tab. B2. 
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Table Bl. Parameter values to be used in eq.s B1-B3 
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Table B2. Parameter values to be used in eq.s B7-B8 
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